Multistability and Critical Fluctuation in a Split Bose-Einstein Condenstate 
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By using a two-mode description, we show that there exist the multistability, phase transition and 
associated critical fluctuations in the macroscopic tunneling process between the halves of a double- 
well trap containing a Bose-Einstein condenstate. The phase transition that two of the triple stable 
states and a unstable state merge into one stable state or a reverse process takes place whenever 
the ratio of the mean field energy per particle to the tunneling energy goes across a critical value 
of order one. The critical fluctuation phenomenon corresponds to squeezed states for the phase 
difference between the two wells accompanying with large fluctuations of atom numbers. 

PACS number(s): 03.75.Fi, 05.30.Jp, 32.80.Pj 

Over the last few years, there has been considerable 
interest in the subjects related to the macroscopic tun- 
neling processes between the halves of a double- well trap 
containing a Bose-Einstein condenstate (BEC). Here we 
only mention a few ones closely related to the topic con- 
cerned in the present paper, for instance, experimen- 
tally realized squeezed states pj, the scheme for demon- 
strating nonlinear Josephson-type oscillations of a driven, 
two-component BEC [|J , coherent oscillations concerning 
Josephson effects, tt oscillations, and macroscopic quan- 
tum self-trapping [^Q. Javanainen and Ivanov fl have 
claimed that these studies are based on essentially 
classical models. On the other hand, the fluctuations 
of atom numbers and phases have been intensively in- 
vestigated quantum-mechanically mainly by two seem- 
ingly quite different approaches: two-mode approxima- 
tion J5|-Q| and the one based on taking atom number and 
phase difference as conjugate quantum variables [p|jTo|. 
The corresponding investigations in an array of traps con- 
taining BEC have also been carried out recently 
Although extensively studied, there still exist some im- 
portant open issues in the the macroscopic tunneling pro- 
cesses between the halves of a double- well trap containing 
a BEC. Some of these open issues are 1) the fluctuations 
of atom number and phase difference obtained by the 
two above mentioned approaches seem to show large dif- 
ference [p 11 1; 2)there exists no any bridge to connect 
the strong and weak tunneling regimes yet JO]; 3) how 
to relate the essentially classical models [JM and the 
corresponding quantum-mechanical description in deal- 
ing with the fluctuations. In this paper, we shall solve 
these three important issues by providing a united ap- 
proach. What is more important, we shall show that the 
previous studies have missed the phenomena of the multi- 
stability, phase transition and associated critical fluctua- 
tions in the macroscopic tunneling processes between the 
halves of a double- well trap containing a BEC. The criti- 
cal fluctuation phenomenon corresponds to the squeezed 



states for the phase difference between the two wells with 
extremely large fluctuation of atom numbers. 

We consider a model system of the many-atom ground 
state of a Bose-Einstein condensate in a double-well po- 
tential p|~|7|,p^|-^6[ . In this system, TV bosonic atoms are 
confined by a infinite harmonic potential that is divided 
into left and right wells by a barrier that can be raised 
and lowered arbitrarily. Making a simple two-mode ap- 
proximation, considering only the lowest energy states, 
the creation and annihilation operators (a L R and clt,.r, 
respectively) for atoms localized in the ground state of 
either the left or the right potential well can be con- 
structed. Neglecting terms that depend only on the total 
conserved particle number N, the Hamiltonian for the 
system can be written [0 



H = — {n L ~ n R ) + —{n 



+ "a) + l{a\a R + a R a L ), 



(1) 



where ht,R = a L R aL,R> the total conserved particle 
number operator N = ul + g = ATra sc h 2 / m is the 
mean-field enegy constant (a sc is the s-wave scattering 
length). The term in 7 describes tunneling between wells, 
whereas the term in g(3, which depends on the number 
of atoms within each wells, describes mean field energy 
due to interactions between atoms in the same well. The 
term in A = {El — Er)/Ti describes the energy differ- 
ence of the ground states in the left and right wells. The 
coefficients 7 and (3 are determined from integrals over 
single-particle wave functions jl],[|. 

Introducing phase operators 4>l,r, by the relation 



[§ exp{4 L , R ) = (n L , R + l)- l l 2 a L ,R or 



aL,R = exp(i<f> L . R ) ^/h L ,R 



(2) 



where we have made use of the fact aL,RF(nL,R) = 
F(fiL,R + ^)o-l,r- The phase and atom number opera- 
tors satisfy the commutative relations, 
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nL,R, <PL,R 



(3) 



But phase operators thus introduced suffer from 
the well-known non-Hermitian problem |l7]-^] that 
(exp(«^ Lj _ R )) t exp(i0 Ljfl ) = 1 - |0)(0| ^ 1 although 
exp(i(j>Lji)(exp(i(j)i J f!)) ji = 1. However, close inspection 
of this problem, we realize that this problem is in fact 
avoidable if we focus on system's states with ni, r =/= 0, 
which will be assumed to be so hereafter. Let n = til and 
<f> — 4>l — 4>r denote the left well's atom number operator 
and the phase difference between the left and right wells 
respectively. Utilizing hr — N — h, and the total atom 
number N = N (the conserved operator N only takes a 
unique eigenvalue N and hence we need not to consider 
its operator characteristic), we can, after omitting the 
unimportant conserved quantity N(A + g/3N)/2, rewrite 
the Hamiltonian (1) as follows, 

H = An — g/3n(N - n) + 7 (\/n(N - h) exp(i</>) + h.c.^j 

(4) 

where the left well's atom number operator n and the 
phase-difference operator tf> satisfy the commutative re- 
lation 



n, <p = 1 



(5) 



Equations (4) and (5) are the fully quantum-mechanical 
model give the united description for investigating the 
various problems in a Bose-Einstein condensate in a 
double-well potential, particularly those relevant to the 
phase coherence. In particular, we shall demonstrate that 
all the previous models, either essentially classical or the 
fully quantum-mechanical ones, dealing with this systems 
can be derived from these equations under some approx- 
imations. 

Suppose that the system can be described by a two- 
mode coherent state \^) = |az,,afl) c characterized by 
two complex parameters oll — y/nexp(i4>L) an d aR = 
\/ N — nexp(i0jj), we can derive from (4) Hamiltonian 
formalism for the equations of motion for the mean atom 
number n— (^\n\^) in the left well (as well as the mean 
atom number N—n in the right well) and phase difference 
4> = ('I' 1 0| 'J) between the two wells as follows (obtained 
by taking average operations to the Heisenberg equation 
of motion dA/dt = i[H, A] for A — n and <j> respectively), 



dn _dU d<(> _ dU 
dt dej) ' dt dn 



(6) 



U = H + An - g(3n(N -n) + 2j v / n(N - n) cos cj> (7) 

where H = and H = N(A + gf3N)/2 is con- 

served quantity and can be omitted without loss of gen- 
erality. The explicit form of equation(6) reads 



^ = 2j v / n(N ~ n) sin </> (8a) 



^ = A-g{3(N-2n) +1 N 2n cos* (8b) 



These are nonlinear versions of the usual Josephson- 
junction equations and are nearly identical to those 
describing the double- well tunneling problem in the same 
system H|]]. Raghavan et al. have investigated Joseph- 
son effects, 7r oscillations and macroscopic quantum self- 
trapping for this system However, no one, to the 
best of our knowledge, seems to have so far noticed the 
important phenomena of multistability and phase tran- 
sition in the macroscopic tunneling process between the 
halves of a double-well trap containing a Bose-Einstein 
condenstate, which we now turn to investigate. 

The multistability and phase transition in the macro- 
scopic tunneling process are clearly seen from Fig. 1. Let 
us describe their main features. First of all, there are two 
kinds of evolution pattern for atom number n (N — n) 
in the left (right) well and phase difference between the 
two wells. They are: 1) stable and unstable steady states 
(fixed points in the phase diagrams) denoting there exists 
no tunneling at all although tunneling rate is non-zero. 
2) periodic tunneling processes where the tunneling am- 
plitude can be very large, i.e., macroscopic quantum tun- 
neling, even for small tunneling rate 7. Secondly, the 
phase transition takes place when an "order" parameter 
|£| = \g(3N/ (2^)1 characterizing the relative magnitude 
of the mean field energy per particle and the tunneling 
energy goes across the critical parameter £ c of order one. 
In other words, one of the three stable fixed point re- 
mains while the unique unstable fixed point and the two 
of the three stable fixed points for |£| > £ c merge into one 
stable fixed point when |£| goes from below to above the 
critical value £ c . Therefore the phase transition corre- 
sponding to the sudden structural change in n — <f> phase 
diagrams when the "order" parameter |£| goes across its 
critical value £ c . Thirdly, There exist three stable and one 
unstable fixed points when |£| > £ c , whereas there exist 
two stable fixed points and no unstable one otherwise. 
The concrete value of the critical parameter depends on 
the parameter 6 = A/(2j) corresponding to the ratio of 
the ground energy difference of the two wells to the tun- 
neling rate. When |£| > £ c , the two of the three stable 
fixed points (i.e., the two stable fixed points at <f> = in 
the lower phase diagrams of Fig.l) are symmetric about 
n = N/2 if the right and left wells are identical to each 
other (5 — 0) and they are asymmetric if the two wells 
are different from each other (8 7^ 0). 
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FIG. 1. The 7Y/(27) contours for several values of param- 
eters £ = gf3N/(2-y) and S = A/ (27). The horizontal and 
vertical axes denote reduced atom number x — n/N and the 
phase difference 4> respectively. The Parameters for the four 
diagrams are (al)(£ = 0.6,5 = 0), (a2) (£ = 1.8,8 = 0), (bl) 
(£ = 1.1, S = 0.1) and (b2) (£ = 1.8, 8 = 0.1) respectively. 

The fixed points and the critical parameter £ c are eas- 
ily shown to be determined by the equations sin </> = 
or cf>Q — 0, 7r and 

(1 - 2x ) (cos</>o ~ 2£- v /x (l ~ ^o)) = 6y/x (l - x ), 

(9) 

where £ = gfiN/{2i), 6 = A/(2 7 ), and x = n /N. 
We find from this equation that the critical parameter 
£ c = 1 for <5 = 0, and £ c > 1 for nonzero 6. For in- 
stance, £ c « 1.37 for 5 = 0.1. In the case 6 = denoting 
identical ground energy for two wells, we can easily ob- 
tain the explicit expressions for the fixed points. They 
are the fixed point P whose no = N/2 and cj>Q = for 
£ > and 0o = tt for £ < 0, the fixed points S± whose 
n = AT(1 ± y / l-£- 2 )/2 and O = for £ > and 
4>o = tt for £ < 0, the fixed point S whose no = N/2 and 
4>o = 7T for £ > and (fo = for £ < 0. The fixed point P 
is unstable for |£| > 1 and stable for |£| < 1, 5 is stable 
and exists for any £, and S± are stable but they exist 
only when |£| > 1. 

Now let us illustrate that (4) and (5) provide a natu- 
ral basis for describing quantum-mechanically the atom 
number and phase statistics. To this goal, we expand 
the left well's atom number operator n and the phase- 
difference operator (f> in the Hamiltonian (4) around one 
of the stable steady states discussed in the paragraph 



where (9) locates, i.e., <fi — no + fj and ij> = <f>o + ip 
with (no,</>o) denoting one of the stable steady states 
(note sin</>o = 0) and [fj, ip] — i. Then after ne- 
glecting the terms equal to or higher than the order of 

O (fj 3 , fjip 2 , fj 2 %p, '0 3 ^ and taking ip — —id/di] and fj = n 
in order to satisfy [fj, ip] = h we can write (4) as follows, 



d_ 

drj 



d_ 

'dr, 



Ejd 2 
2 drf 



E iiz + E wi--^-i-o + ^v 2 (10) 



where Hq is a constant equal to (H — Ho) in (7) evalu- 
ated at the fixed point considered, and the coefficients in 
(10) are given in (11) for the fixed point S and in (13) 
for the fixed points S± respectively. Following the same 
argument as the one in ref. j^] given in the paragraph im- 
mediately after its equation (16) , we can neglected £1,2 
terms in (10) in evaluated the atom number and phase 
statistics. After this approximation, (10) is nothing but 
a harmonic oscillator model and fluctuations of the atom 
number in the left well and the phase difference are eas- 
ily shown to be An = (A0)" 1 = (Ej/E c ) 1/4 when the 
harmonic oscillator is in its ground state. 

In the case where 6 = denoting identical ground en- 
ergy for the two wells, if one considers the stable fixed 
point S whose concrete expression is given after (9), the 
coefficients in (10) can be calculated to be 



F F gfj F 9l3N2 



En = 



2g/3 



(11) 



Consequently, fluctuations of the atom number in the left 
well and the phase difference in this case are given by 



(12) 



a 1 VN 

A7^ = -7— = — = 

A0 V2(l 



where £ = g(3N/(2~f). The equation (16) of ref. g deal- 
ing with the same system as ours has nearly identical 
form as our (10) and (11) but with two slightly different 
coefficients E c = 2g(3(l + 0.5 |CI ~ 1 ) and E 2 = 0.5g/3/|£| 
in our notation. It is instructive to note that the atom 
number fluctuation An oc iV 1 / 4 in the weak tunneling 
regime (|£| ^> 1) just as given by Leggett and Sols M, 
whereas in the strong tunneling regime (|£| <C 1), the 
atom number fluctuation An oc y/W just as obtained by 
Javanainen and Wilkins [0,0- Therefore the problem 
quarreled by them |8|,[ll| are naturally settled. 

The previous studies on the atom number statistics for 
the double- well system have failed to notice the multista- 
bility and the phase transition in this system when the 
parameter |£| goes across the critical parameter £ c . Con- 
sequently, no one has so far discussed the atom number 
and phase statistics around the stable fixed points S± 
when |£| > £ c , which will be the subject of the present 
paragraph. The critical parameter £ c is unity in the case 
where 6 = denoting identical ground energy for the two 
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wells. In this case, we consider fluctuations around the 
stable fixed points S± whose concrete expression is given 
after (9), and obtain the coefficients in (10) as follows, 

E x = ±0.5 5/ 3iVyi-£- 2 , E 2 = -gfc 2 , (13a) 
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E c = -2g/3(e - 1). 



(13b) 



Fluctuations of the atom number in the left well and the 
phase difference in this case are given by 



An 



Acj) 



(14) 



where |£| = \g(3N/ (2j)\ > 1. In the weak tunneling 
regime (|£| 1), the atom number fluctuation Sn ss 
is very small in large-N circumstances, which 
demonstrates sub-Poissonian fluctuations and the atoms 
in any one of the wells can be thought to be approxi- 
mately in a Fock state. Another interesting phenomenon 
is that the atom numbers in both wells display strong 
fluctuations when the parameter |£| = \gf3N/(2j)\ ap- 
proaches one from above. The atom number fluctuations 
have the form An w y/N/2[2(\£\ - 1)]~ 1/4 as |£| -> 1. 
This form demonstrates the typical strong critical fluctu- 
ation phenomena in phase transitions [ pT|]2^ ] . The corre- 
sponding "order" parameter and critical index for atom 
number fluctuation in our case are |£| and 1/4 respec- 
tively |2^j2^]. However, It should be emphasized that the 
critical fluctuation phenomenon corresponds in fact to 
squeezed states for the phase difference cj> since the phase 
difference fluctuation A0 w y/2/N[2{\£\ - l)] 1 / 4 -> as 
the "order" parameter £ approaches unity from above. 

The multistability and phase transition as well as 
the critical fluctuation phenomenon (i.e., phase squeezed 
states) in the macroscopic tunneling in a double- well trap 
containing a Bose-Einstein condenstate(BEC) are well 
within the reach of nowadays BEC-related technology. 
As a matter of fact, the recent experiment by Kasevich's 
group JjJ has already reached the strong tunneling regime 
and the parameter £ = gf3N / (27) in that experiment can 
at least reach as low as 1.5 as given in the caption of 
its Fig. 1 D 0. We therefore believe that one should 
be able to discover the new phenomena investigated here 
with the same apparatus as in that experiment. 
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